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1. INTRODUCTION 

■ 

o ■ 

Coulomb gauge Yang-Mills theory (and by extension Quantum Chromodynamics) is a fascinating, yet frustrating 
Cn| . endeavor. On the one hand, Coulomb gauge offers great potential for understanding such issues as confinement [1, 2]; 
on the other, the intrinsic noncovariance of the formalism makes any perturbative calculation formidably complicated. 
Many approaches to solving (or providing reliable approximations to solving) the problems in Coulomb gauge have 
been forwarded. Recent among these are the Hamiltonian approach of Ref. [3] , based on the original work of Christ 
and Lee [4] . A lattice version of the Coulomb gauge action also exists [5] , which has led to numerical studies, for 
i example Refs. [6]. Functional methods based on the Lagrangian formalism have also been considered, especially within 
the first order (phase space) formalism [1, 7] and most recently, one-loop perturbative results for both the ultraviolet 
divergent and finite parts of the various two-point functions have been obtained [8]. Similar results were previously 
obtained for the gluon propagator functions under a different formalism (using the chromoelectric field directly as a 
degree of freedom and without ghosts) and using different methods to evaluate the integrals [9] . 

In this paper, we consider the (standard, second order) functional approach to Coulomb gauge Yang-Mills theory. 
We derive the Dyson-Schwinger equations and Slavnov-Taylor identities for the two-point functions that arise in the 
construction and using the techniques of [8] we present results for the one-loop perturbative dressing functions. 

The paper is organized as follows. In the next section, the functional formalism used is described. Section 3 
concerns the decomposition of the functions used. The (nonperturbative) Dyson-Schwinger equations and Slavnov- 
Taylor identities relating the various Green's functions are derived in Section 4. In Section 5, the one-loop perturbative 
\Q , results are obtained. Finally, there is a summary and outlook. 

o\ 
m 

! 2. FUNCTIONAL FORMALISM 

o 

Let us begin by considering Coulomb gauge Yang-Mills theory. We use the framework of functional methods 
to derive the basic equations that will later give rise to the Dyson-Schwinger equations, Slavnov-Taylor identities, 
Feynman rules etc. Throughout this work, we will use the notation and conventions established in [7, 8]. We work 
k^j ' in Minkowski space (until the perturbative integrals are to be explicitly evaluated) with metric = diag(l, — 1). 
Greek letters (fi, v, . . .) denote Lorentz indices, roman subscripts (i, j, . . .) denote spatial indices and superscripts 
(a, b, . . .) denote color indices. We will sometimes also write configuration space coordinates (x, y, . . .) as subscripts 
where no confusion arises. 

The Yang-Mills action is defined as 



Sym = I d 4 x 



(2.1) 

where the (antisymmetric) field strength tensor F is given in terms of the gauge field A a : 



"4 H v 



F% = d » A ° - d »K + 9f abc A b a A c u . (2.2) 

In the above, the f abc are the structure constants of the SU(N C ) group whose generators obey [T a ,T b ] = if abc T c . 
The Yang-Mills action is invariant under a local SU{N C ) gauge transform characterized by the parameter 9%: 

U x = exp{-i9 a x T a }. (2.3) 

The field strength tensor can be expressed in terms of the chromoelectric and chromomagnetic fields (a = A ) 
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E a = -d a A a - Va a + <?/ ahc iV, Bf = e ijk 



VjA a k ^9.f abc A)At 



(2.4) 
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such that Sym = J {E 2 — B 2 )/2. The electric and magnetic terms in the action do not mix under the gauge transform 
which for the gauge fields is written 



A» -+ 4 = U X A,U X - - g {d^U x )Ul 

Given an infinitesimal transform U x = 1 — iQ%T a , the variation of the gauge field is 

1 



_ D acQC 



(2.5) 



(2.6) 



(2.7) 



where the covariant derivative in the adjoint representation is given by 

Dl c = 5 ac d^ + gf abc Al. 
Consider the functional integral 

Z = J V<f>cxp{iS Y M} (2.8) 

where $ denotes the collection of all fields. Since the action is invariant under gauge transformations, Z is divergent 
by virtue of the integration over the gauge group. To overcome this problem we use the Faddeev-Popov technique 
and introduce a gauge- fixing term along with an associated ghost term [10]. Using a Lagrange multiplier field to 
implement the gauge-fixing, in Coulomb gauge (V • A = 0) we can then write 



£>$ exp {iSym + iS f P }, S 



fp 



d 4 x 



-X a V-A a - c a V-D ab c b 



(2.9) 



The new term in the action is invariant under the standard BRS transform whereby the infinitesimal gauge parameter 
9 a is factorized into two Grassmann-valued components 9 a = c a 5X where 5X is the infinitesimal variation (not to be 
confused with the colored Lagrange multiplier field A a ). The BRS transform of the new fields reads 



Sc a 

8c a 
8X a 



-X a SX 
9 

_l f abc c b c c sx 

2 

0. 



(2.10) 



It is at this point that this work diverges from Rcf. [7] in that we remain here within the standard (second order) 
formalism. By including source terms to Z, we construct the generating functional, Z[J]: 



Z[J] = J V<S> exp {iSym + tSfp + iS s } 



where 



/ 



rf 4 



p a a a + J a -A a + cV + rj a c a + £ a X a 



(2-11) 



(2.12) 



It is convenient to introduce a compact notation for the sources and fields and we denote a generic field $ Q with 
source J a such that the index a stands for all attributes of the field in question (including its type) such that we can 
write 



J a <S> a 



(2-13) 



where summation over all discrete indices and integration over all continuous arguments is implicitly understood. 
Expanding the various terms we have explicitly 



Sym - j d A x j-^f [5 l3 d 2 - «%V 2 + VjV,-] A{ - A{d V i( r f - X -a^ 2 a^ 
+9f fbc \-(d A{)A b a c - {V i( jf)A\cj c + (VjAl^At] + g 2 f fbc f fde 



\A\a c A d X - \A b A)AtA] 



(2-14) 
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The field equations of motion are derived from the observation that the integral of a total derivative vanishes, up 
to boundary terms. The boundary terms vanish, although this is not trivial in the light of the Gribov problem [2] 
(the reader is directed to Ref. [7] and references therein for a discussion of this topic). Writing S — Sym + Sf P , we 
have that 







V1>—-exp{iS + iS s }. 



(2.15) 



The explicit form of the field equations of motion is given in Appendix A. 

In addition to the field equations of motion, there exist identities derived by considering the BRS invariance of the 
action (these eventually form the Slavnov-Taylor identities). The BRS transform is continuous and we can regard it 
as a change of variables in the functional integral. Given that the Jacobian of such a change of variables is trivial and 
that the action is invariant, we have that 
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(2.16) 



= / — exp {iS + iS s + iSS s } sx=0 



J X>$exp{tS + tS a } J d 4 



-p a d c a + f abc p a a b c c - -jy iC a + f abc J?A b c c + -\ a rf + -f abc r] a c b c c 
9 9 9 2 



So far, the generating functional, Z[J], generates all Green's functions, connected and disconnected. The generating 
functional of connected Green's functions is W[J] where 



Z[J] 



,W[J] 



We define the classical fields to be 



4\ 



hi 



X>$ $ Q exp iS 



1 SZ 
Z 8iJ a 



(2.17) 



(2.18) 



The generating functional of proper Green's functions is the effective action, T, which is a function of the classical 
fields and is defined through a Legendre transform of W: 



r[$] = w[j) - ij a <, 



(2.19) 



We introduce a bracket notation for derivatives of W with respect to sources and of V with respect to classical fields 
(no confusion arises since the two sets of derivatives are never mixed) : 



<iJ a >- 



SW 
5iJ a 



<l$ a >-- 



ST 



(2.20) 



It is now possible to present the field equations of motion in terms of proper functions (the Dyson-Schwinger equations 
are functional derivatives of these equations). Using the results listed in Appendix A we have: 

<iA a ix > = - [5ijdl x - SijVl + V ix V jx ] A% - d ax V lx a a x + V ix X a x 

+9f abc J d A yd A zd 0x 5(y - x)S(z - x) [<iJ b y ip> +A b y a c z ] 

~9f fac J d A yd A z5{z - x)V ix S(y - x) [<ip{ip c z > + <vrf z tr]l> +a(a c z + c f y c z ] 

+9.f abc j ' d A yd A z [SijS(z - x)\7 kx S(y - x) +S jk 6(y - x)V ' ix 5(z - x) - S ki V jx S(y - x)S(z - x)] x 

[<ljj y lj kz > +A b y Al z ] 

+g 2 f fac f fde [ <lp c xl jf x ^> +(J c <ljf xl p x> +<J X <ip x lJ? x > +Af x <l Px l Px > +<7^iX] 



-lg2ffbcffde 6jk6u [ 6 gc § eh( § ab § di + gad^ + gbg gdh (gacgei + gaegci^ x 



[<lJ 9 jx ljtlJl> +A% <%J^tj} x > +A\ X <lJ 9 rx lj£ x > +4* <*Jjx*4c> + A % A kxMx] 

<ia a x > = -d Qx V ix A a ix - V 2 x a a x - gff ba J d 4 y d 4 z S(z - x)d 0x S(y - x) 



(2.21) 



<lj( y lJ iz > +A{ y A iz 
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+9f abc J d 4 yd 4 z [V ix 8{y - x)S(z - x) + 5(y - x)V ix 5(z - x)} [<iJ h iy ip> +A b y a c z ] 

+9 2 f fba f fde [<iJ h ix iJt x ip> +A\ X <iJ? x ip e x > +a e x <iJ b x iJf x > +Af x <%J b ix ip x > +A b ix Ala c x 
<iX a x > = -V ix A a ix , 

<«■> = -V 2 x c a x + gf abc J d 4 yd 4 zV> x 6(y - x)5(z - x) [<iJ b y iff z > +A b y 4] ■ 
It is also useful to express the A equation of motion in terms of connected functions: 

C = v ix <iJ? x > . 



(2.22) 
(2.23) 

(2.24) 



(2.25) 



The identity stemming from the BRS invariance is also best expressed in terms of both connected and proper functions 
and reads: 



= J d 4 x <t Q> +- gP a x 8 0x <irf x > +f abc P : [<ip h x irf x > + <ip x Xirf x >\ - ± 

w + <<><«>] + [<**> + <*><*>] } • 

= J d^[-~ g <<> K-~ g <<> d 0x c a x - f abc <%a a x > [<ip b x irf x > +^ b x c x ] - - g 

-/* <,a L > u, w [<^«> + J/* «> [<•*«> +44 } , 



»z ja 



(-v|) 



(2.26) 



(-V2) 



<tc"> 



(2.27) 



where we have used the common trick of using the ghost equation of motion in order to reexpress one of the interaction 
terms transversely, with the transverse projector in configuration space being Uj(x) = Sij + ^7i X ^7j X /{— V^). This 
manipulation will be useful when we consider the Slavnov- Taylor identities for the two-point functions later on. 

At this stage it is useful to explore some consequences of the above equations that lead to exact statements about 
the Green's functions. Introducing our conventions and notation for the Fourier transform, we have for a general 
two-point function (connected or proper) which obeys translational invariance: 



<iJ a (y)iJfj(x)> = <iJ a (y — x)iJp(<S)>= j dkW a p{k)e~ 
<i$> a (y)i$p(x)> = <i$ a (y-x)i$p(0)>= J dkT al3 (k)e 



ik-(y-x) 

i 



where dk — d 4 k / (2ir) 4 . Starting with Eq. (2.23), we have that the only non-zero functional derivative is 



<iA b y i\ a x >= i5 ba V ]X 8{y -x) = 5 ba J dkk je - %k < y -^ 



(2.28) 



(2.29) 



and all other proper Green's functions involving derivatives with respect to the A-ficld vanish (even in the presence 
of sources). In terms of connected Green's functions, Eq. (2.23) becomes Eq. (2.25) and the only non-zero functional 
derivative is 



V fa <i£iJ? x >=-i5 ba 6{y-x). 



(2.30) 



Because Eq. (2.25) involves the contraction of a vector quantity, the information is less restricted than previously. 
However, we can write down the following (true once sources have been set to zero such that the tensor structure is 
determined): 



<i3 b v i3 a lx > = J dkW b A a A (k)U,(k)e 



-ik-(y—x) 



<%p b y lJt x > = 0, 



(2.31) 



where tji(k) 



kjki/k 2 is the transverse projector in momentum space. These relations encode the transverse 



nature of the vector gluon field. Turning to Eq. (2.26), we recognize that if we functionally differentiate with respect 
to iijy, again with respect to i£| and set sources to zero, we get that 



(2.32) 



■5 



w m 




a 


A 


Ai 







(-hi) 


tij (k) 2 r*2\ 
' _ ' 


fc2 


a 





p 


(-fc°) 


A 


kj 

p 


p 


lo 



pv u ,< 


Aj 


a 


A 


A, 


iktyS-ij ~\~ %k tijtk*} 


%k k^i 


k't 




%k kj 


-ik 2 





A 


—kj 









TABLE I: Tree-level propagators [left] and two-point proper functions [right] (without color factors) in momentum space. 
Underlined entries denote exact results. 

In effect, the auxiliary Lagrange multiplier field A drops out of the formalism to be replaced by the transversality 
conditions, as it is supposed to. 



FEYNMAN RULES AND DECOMPOSITIONS 



Let us now discuss the Feynman rules and general decompositions of Green's functions that will be relevant to 
this work. The Feynman rules for the propagators can be derived from the field equations of motion (written in 
Appendix A) by neglecting the interaction terms and functionally differentiating. Denoting the tree-level quantities 
with a superscript (0), the corresponding equations read: 



J? x = Mo, - SijVl + V ix V jx ] <3Jj> (0) +d» x V lx <% P %><® -V fa <*£>< >, 

P a x = d 0x v iX <iJ? x > {n) +v 2 < v s><°\ 

= V 2 <irf x > (0) ■ 



T a 7n b > (0) = 



The tree-level ghost propagator is then 

<W x m" y 

and we identify the momentum space propagator as 



l S ab [ dk^e- 1 
J fc 2 



fe-(j/-x) 



-5 ab ±. 
k 2 



(3.1) 



(3-2) 



(3.3) 



The rest of the propagators follow a similar pattern and their momentum space forms (without the common color 
factor 5 ab ) are given in Table I. Note that it is understood that the denominator factors involving both temporal and 
spatial components implicitly carry the relevant Feynman prescription, i.e., 



1 



(fc 2 - fc 2 ) (fc 2 - fc 2 + 20+) 



(3.4) 



such that the integration over the temporal component can be analytically continued to Euclidean space. It is also 
useful to repeat this analysis for the proper two-point functions and using the tree-level components of Eqs. (2.21), 
(2.22) and (2.24) we have 



<iA a ix >^ = -[Sijd^-SijVl + V^Vj^A^-do^i^ + V^Xl, 
<ia a x >^ = -do.V^-V 2 ^, 



<zc> (0) = -V 2 x c a x . 



The ghost proper two-point function in momentum space is 

r(°) afc (fc) = 5 ab ik 2 



(3.5) 



(3.6) 



and the rest are presented (without color factors) in Table I. It is immediately apparent that the gluon polarization 
is not transverse in contrast to Landau gauge. 
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w 




a 


A 


Ai 







fc 2 


a 









A 


kj 








r 




a 


A 


a, 


-i{kl - k 2 )tij(k)T aa - ikl^-V A A 


ik kiT^ Act 


k% 




%k kjY'Ao 


lk TcTfj 





A 


—kj 









TABLE II: General form of propagators [left] and two-point proper functions [right] (without color factors) in momentum space. 
All dressing functions are functions of fcp an d k 2 . 

The tree-level vertices are determined by taking the various interaction terms of Eqs. (2.21-2.24) and functionally 
differentiating. Since, in this study, we are interested only in the eventual one-loop perturbative results we omit the 
tree- level four-point functions (1^4 and Taa<jo)- Defining all momenta as incoming, we have: 

r { °l%( Pa ,p b , Pc ) = i 9 r bc 5 jk (p b -p% 



^(Q)abc 
crAcrj 



(Pa,Pb,Pc) 



r( 3Mjk(Pa,Pb,Pc) 

p(0)a&c 
ccAi 



-igf ahc {Pa-Pc) 3 i 

-igf abC [Sij{p a - Pb)k + Sjk(pb - Pc)i + Ski(p c - Pa)]] , 



(Pc,Pc,pa) = -wf a C Pc>, 



(3.7) 



In addition to the tree-level expressions for the various two-point functions (connected and proper) it is necessary 
to consider their general nonperturbative structures. These structures are determined by considering the properties of 
the fields under the discrete transforms of time-reversal and parity (the noncovariant analogue of Lorcntz invariance 
arguments for covariant gauges). Using the same techniques as in Ref. [7] we can easily write down the results in 
momentum space. For the ghost, we have 



Wf{k) = -6 ab —D c {k 2 ), rf (fc) = 5 ab ik 2 T c (k 2 ) 
k^ 



(3.8) 



and the rest are presented in Table II. With the exception of the ghost, all dressing functions are scalar functions 
of two independent variables, k 2 and k 2 . The ghost dressing functions are functions of k 2 only for exactly the same 
reasons as in the first order formalism [7]. At tree-level, all dressing functions are unity. 

The dressing functions for the propagators and two-point proper functions are related via the Legendre transform. 
The connection follows from 



5iJp 
8iJ„ 



-l-J^- <l$f3> = J^r- <l$~ / l$fj> = <lJ a lJ~ f Xl$~ f l$fj> 



8iJ n 



(3.9) 



(Recall here that there is an implicit summation over all discrete indices and integration over continuous variables 
labeled by 7.) Considering all the possibilities in turn, we find that 

d A a = r^ 1 , 
n — r -1 

-1 



D r = r 



c ' 



(3.10) 



Actually, whilst we have included D a \ up to this point, since there is no vertex involving the A-field this propagator 
will not directly play any role in the formalism. However, indirectly it does turn out to have a meaning as will be 
shown in the next section. 



4. DYSON-SCHWINGER EQUATIONS AND SLAVNOV-TAYLOR IDENTITIES 

With the observation that 

fi 

<iJ 7 iJ a >~ — <^J 7 ^Je><^$ £ ^$/3^$^><^J^^J a > 



5i$ 



(4.1) 
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[stemming from the Legendre transform and following from Eq. (3.9)], the derivation of the Dyson-Schwinger equations 
becomes relatively straightforward. Starting with Eq. (2.21), omitting the terms that will not contribute at one-loop 
perturbatively and recognizing the tree-level vertices in configuration space, we have that 

<iA a ix > = ilSijd^-SijVl + VixVj^iA^+idoxVixiaZ-iVixiXl 

-J d 4 yd 4 zT^ A %(z,x,y) [<tJ b y ip c z > -%A) y ujf[ -J d*yd 4 z ^T^(z,x,y) [<ip b y zp> -ia b y ia<;] 

-J d*ydH±T^ c k {x,y,z) [<iJ b y zJ c kz > -^A c kz ] + J d 4 yd 4 zT^ a (y,z, X ) [<vff z r, b y > +ic* z ic b y \ 
+ ... (4.2) 

Taking the functional derivative with respect to iAf w , using Eq. (4.1), setting sources to zero and Fourier transforming 
to momentum space (each step is straightforward so we omit the details for clarity) we get the Dyson-Schwinger 
equation for the gluon polarization: 

T AAii( k ) = Saf [~< k t - k 2 )S t i - ikiki 

+ J { J tfwrStSV-fc,A.-w)W^Hlf^«(w,-fe,fc-a;)^(a;-fc) 

+i J ^ r£(i, w - fcX^j^rf^jc, -fc, fc - - fc) 

- / lS> - k, -co, k)W* d (u)r*? Al (u, k co, -k)Wf{co -k) + ... (4.3) 

Turning now to Eq. (2.22), we have 

<<> - id 0x V ix iAl + iV 2 x ia a x - J d 4 yd 4 z ^vf A abc jk {x,y, z) [<iJ b y iJ c kz > -iA b y iA c kz ] 

- J d 4 y d 4 z T^{x, y, z) [<iJ b jy ip c z > -iA b y ia z ] +... (4.4) 

where again, terms that do not contribute at the one-loop perturbative level are omitted. There arc two functional 
derivatives of interest, those with respect to zct£ and %A[ w , which give rise to the following two Dyson-Schwinger 
equations: 

Ki(k) = + I J a W r^ fc fo-a;,a;-*0^ 

+ / ^r<°]J(fc, - W , W - i=X Aim ( W )r^ m (-fc, W , - - fc) + . . . (4.5) 

+ J dcoTf A ab ;{k,-co,co-k)W b A 3 MT^ (4.6) 
Next we consider the ghost equation, Eq. (2.24), which can be written 

<«S>= iV*u£ + J d 4 y d 4 z T^°(x, y, z) [<iJl z irf> y > -iA\ z ic\\ . (4.7) 
The ghost Dyson-Schwinger equation is subsequently 

T a J{k) = 5 a hP + J iwr^(k,-u>,w-k)Wy(u)r^.(u,-k,k-u)W (4.8) 



In addition to the Dyson-Schwinger equations, the Green's functions are constrained by Slavnov-Taylor identities. 
These are the functional derivatives of Eq. (2.27). Since Eq. (2.27) is Grassmann- valued, we must first functionally 
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differentiate with respect to ic d . We are not interested (here) in further ghost correlations, so we can then set ghost 
sources to zero. Also, there is no further information to be gained by considering the Lagrange multiplier field A a , 
and we set its source to zero also. Equation (2.27) then becomes 



-d Qy <ia d > -f abd <ia a y > ia b - f abd iA b y t 3l {y) <iA a iy > 

= |^|_ / a fcc<<> _^ < ^- >+ I 



(-v|) 



<iA?> 



<ic a x ic d > -f abc <iA a ix > Uj{x)-^ <iJ b x vq c x >^ . 

(4.9) 



Taking the functional derivatives of this with respect to io e z or iA% z and setting all remaining sources to zero gives 
rise to the following two equations: 



-d ay <ia e z ia d > 



-d 0y <iA e kz ia d > 



v ix e , a . 

(_V 2 ) Z lx 



<ic a x ic d > 



/*•{; 

-r hc <icT e z ia a x > <ip h x vr&> -f abc <iatiA1 x > U 3 (x)-^ <oJ^>} , (4.10) 



d 4 x i - 



l kz 1 ix" 



<ic x ic y > 



-f abc <iA% z ia a x > <ip h x irf x > -f abc <iA% z iA a ix > k 3 {x)-^ <iJ b x vff x >^ 



SlCy 



(4.11) 



Now, using Eq. (4.1), we have that 



f abC J^d <ip b x vrfx>= -f abc <i^ x ir] a ><ic a ic d y i^ 1 ><iJ 1 ip x >= 5 ad t a - cc {x,y). 



(4.12) 



Taking the Fourier transform 



^a;cc{ x , U) — J 



d k S CT; c C (fc)e 



-%k-{x— y) 



(4.13) 



we get that 



du> W c (k — uu)T- cl {k - oj, —k, uj)W ia {uj). 



(4.14) 



Since the ghost Green's functions are independent of the ghost line's energy scale [7], after w has been integrated 
out, there is no external energy scale and 



(4.15) 



However, under time-reversal the cr-ficld changes sign (such that the action remains invariant) which in momentum 
space means that under the transform k a — ► — k , S CT; c C (fc) must change sign and so, given Eq. (4.15) we have the 
result that 



^a;cc{k) — 0. 



(4.16) 



In the case of the term 



S a f£ Aj - cc (x,y) = J ahc — d <iJ b ]X irf x > 



(4.17) 



we can see automatically that in momentum space, Y*Aj;- c (k) <~ kj and that the transverse projector that acts on it 
in Eqs. (4.10) and (4.11) will kill the term. We thus have 



-do y <io%ia„> 



d 4 x {- 



-d 0y <iA% z ia d > 



/<*{; 



<iA% z iA a ix > 



<<lCy> 



(4.18) 
(4.19) 



9 



which in terms of the momentum space dressing functions gives 

r CTCT (fc 2 ,fc 2 ) - T Aa (k 2 ,k 2 )T c (k 2 ), (4.20) 

r\4 CT (fc 2 ,fc 2 ) - r^(fc 2 ,fc 2 )r c (fc 2 ). (4.21) 

The Slavnov-Taylor identities for the two-point functions above are rather revealing. They are the Coulomb gauge 
equivalent of the standard covariant gauge result that the longitudinal part of the gluon polarization remains bare 
[11]. We notice that they relate the temporal, longitudinal and ghost degrees of freedom in a manner reminiscent of 
the quartet mechanism in the Kugo-Ojima confinement criterion [12]. Also, they represent Gau8' law as applied to the 
Green's functions. Equation (4.9) suggests that proper functions involving the temporal cr-field can be systematically 
eliminated and replaced by functions involving the vector A and ghost fields although whether this is desirable remains 
to be seen. 

We can now return to the general decompositions of the two-point functions. We see that as a consequence of either 
of the two Slavnov-Taylor identities above, Eqs. (4.20) or (4.21), Eq. (3.10) reduces to D a \ = D c , reassuring us that 
at least the formalism is consistent. We also see that there are only three independent two-point dressing functions, 
whereas (accounting for the tensor structure of the gluon polarization) we have five Dyson-Schwinger equations. We 
will investigate this perturbatively in the next section. 



5. ONE-LOOP PERTURBATION THEORY 



Let us now consider the one-loop perturbative form of the two-point dressing functions that are derived from the 
Dyson-Schwinger equations. So far, all quantities are expressed in Minkowski space. The perturbative integrals must 
however be evaluated in Euclidean space. The analytic continuation to Euclidean space (kg — > ikn) is straightforward 
given the Feynman prescription for denominator factors. Henceforth, all dressing functions will be written in Euclidean 
space and are functions of k\ and k 2 . The Euclidean four momentum squared is k 2 = k\+k 2 . We write the perturbative 
expansion of the two-point dressing functions as follows: 

T aP = l + g 2 T^. (5.1) 

The loop integrals will be dimensionally regularized with the (Euclidean space) integration measure 



du>4 d d oj 



(spatial dimension d — 3 — 2s) . The coupling acquires a dimension: 

S 2 ->ffV, (5-3) 

where /U is the square of some non-vanishing mass scale squared. This factor is included in such that the 
new coupling and are dimensionlcss. By inserting the appropriate tree-level factors into the Dyson-Schwinger 
equations, extracting the color and tensor algebra we get the following integral expressions for the various two-point 
proper dressing functions: 

(rf-i)r«(fc 2 ,fc 2 ) = -n c [ M £ ^(fc4 + ^ 4 ) 2 g f ^du: — 

AA J k 2 LJ 2 (k-L0) 2 J k 2 L0 2 (k-L0) 2 

-2N C J ^2 J^fc _ u y tii(k)tjm{u)tnk(k - Co) [5 i: jkk - SjkUJi - Skikj] [S m ik n - Si n k m - 5 nm ui] , 

(5.4) 



r { AA(kl,k 2 ) = -N c [ M £ ^(fc4 +o; 4 ) 2 f 
AAy 4 ^ J k 2 k 2 u 2 {k-Z) 2 3 3 V ' J k\k 2 Q 2 {k - bjf 



ix e dc 



. 2 



—k-(2u — k) — k-Cok- (Co — k) 

9 Nc I VrTJr ^2 M^)Mfc - (5-5) 

T^(klP) - -\N C J ~_ y tij {Q) tji (k - Co) - 4N C J pf^^ UACo), (5.6) 



1 f fi e duo k-(k — 2lo) 
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/i e dbJ (fc 4 — 2CJ4) 



k4k 2 L*J 2 (k — Uj) 2 

T^{k 2 ) = -A c / ^J£J^L k^UAuj). 

c V ' 7 k 2 u 2 (k-uj) 2 1 A ' 



k-(k- 2u)tij {u)tji {k - u) - 2N C 



-kikjtijiw), (5.7) 
(5.8) 



At this stage, we are in a position to check the two Slavnov-Taylor identities for the two-point functions. The first 
of these, Eq. (4.20), reads at one-loop: 



r (i) _ r (i) _ r (i) = o 



Inserting the integral expressions above and eliminating overall constants, the left-hand side reads 

— 2u))tij{UI)tji(K — U) — I — = — 

k^k 2 uj 2 {k — uj) 2 J k 2 uj 2 {k — uj) 2 



rS - rffi - r<.» ~ -\ / j£^M* - ->M*(*- a) - / 



kikjtij (cj) . 



(5.9) 



(5.10) 



By expanding the transverse projectors and scalar products, it is relatively trivial to show that this does indeed vanish. 
The second identity, Eq. (4.21), reads 



r&-r2i-r« = o 



(5.11) 



and the left-hand side is: 

.(1) _ r (i) 1 f dm k-(k - 2g) 



1 Act ~ 1 — 1 <■ 



w 2 (fc — w) 

L0 2 {k-L0) 2 



2 Hk-totoWtiG-anf 



1-, -. 2 

-fc-(2w — fc) — k-ujk-(oj — k) 



dui (u 2 + 2/C4W4) 

w 2 (fc — (J) 2 



kikjtij (c^) 



(5.12) 



Again, it is straightforward to show that this vanishes. Thus, we have reproduced the Slavnov-Taylor identity results 
that tell us that there are only three independent two-point dressing functions. 

The evaluation of the integrals that give Taa, T aa and T c is far from trivial. However, using the techniques 
developed in [8] it is possible. For brevity, we do not go into the details here and simply quote the results. They are, 
as e —* 0: 



r£>(*,y) 
if>(y) 



A c 



( 47r )2-s 

A c 

( 47r )2- £ 

A c 



1 



7 — In 



/i 



64 



3z + g(z) 



1 14 

Tz ~ y 



/(*) 



i 



- 1 + 112 -3z 2 



}• 



(47T) 



2-e 
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y 

4 
3 



— 7 — In 



7 — In 



31 



-- + 6z + .g(z)(l-3z)-/(z) 



1 „ 3 2 

2 +2z+ 2 Z 



28 8 



In 2 , 



(5.13) 



where x = k 2 , y = fc 2 , z = x/y and we define two functions: 



4 In 2—= arctan \fz — I —= ln(l 

y/z Jo Vt(l + Zt) 



m = 

g{z) = 21n2-ln(l + z) 



zt), 



(5.14) 



(The integral occurring in f(z) can be explicitly evaluated in terms of dilogarithms [8].) Defining a similar notation 
for the perturbative expansion of the propagator functions: 



D a!i = 1 + g 2 D ( ^ 

we then have, via Eq. (3.10), the final results: 

D A 1 A (x,y) = -T^ A (x,y), D%(x,y) = -T^(x,y), D^(y) = -lf)(y). 



(5.15) 



(5.16) 



Several comments are in order here. Firstly, the expressions for and Taa, Eqs. (5.4) and (5.5), respectively, 
contain energy divergent integrals of the form 



/ 



duj {l,U!i,U!jUJj} 

uj 2 (k — uj) 2 



(5.17) 
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These integrals cancel explicitly, though it should be remarked that this cancellation is more obvious in the first order 
formalism [8]. Secondly, with respect to the temporal variable x, all the results above are strictly finite for Euclidean 
and spacelike Minkowski momenta - any singularities occur for z = x/y = —1 (the light-cone) with branch cuts 
extending in the timclikc direction. This means that the analytic continuation between Euclidean and Minkowski 
space can be justified. Thirdly, the coefficient of the e-pole for D aa and the combination DaaD 2 is lliV c /3(47r) 2 
which is minus the value of the first coefficient of the /3-function. This confirms that g 2 D aa [13] and g 2 DAAD 2 (the 
Coulomb gauge analogue of the Landau gauge nonperturbative running coupling) are renormalization group invariants 
at this order in perturbation theory. Fourthly, the results above for Daa, D aa and D c are identical to those calculated 
within the first order formalism [8]. 



The two-point functions (connected and proper) of Coulomb gauge Yang-Mills theory have been considered within 
the standard, second order formalism. Functional methods have been used to derive the relevant Dyson-Schwinger 
equations and Slavnov-Taylor identities. One-loop perturbative results have been presented and the Slavnov-Taylor 
identities that concern them verified. 

Suffice it to say that it is tautological for the situation in Coulomb gauge to be somewhat different from covariant 
gauges such as Landau gauge. The proper A- A two-point function is explicitly not transverse, nor does its longitudinal 
component remain bare beyond tree-level. This longitudinal component can however be written in terms of the 
temporal gluon and ghost two-point functions via the Slavnov-Taylor identities. Indeed, the Slavnov-Taylor identities 
show that there are only three independent two-point dressing functions: the (transverse) spatial gluon propagator 
dressing function (Daa), the temporal gluon propagator dressing function (D CT(T ) and the ghost propagator dressing 
function (D c ). With the exception of the ghost dressing function, all are noncovariantly expressed in terms of two 
variables: fc| (or k 2 in Minkowski space) and k 2 . Perturbatively it is seen that the analytic continuation between 
Euclidean and Minkowski space (and vice versa) is valid and that the Slavnov-Taylor identities hold. 

There are many further questions to be addressed. The perturbative structure of the vertex functions, the addition 
of the quark sector and the construction of physical scattering matrix elements from noncovariant components are 
all important next steps. The issue of noncovariant renormalization prescriptions must also be understood. The 
connection of the functional formalism with other approaches such as the Hamiltonian formalism [3] and lattice 
calculations must also be established. Clearly, there is a lot of work yet to be done. 



This work has been supported by the Deutsche Forschungsgemcinschaft (DFG) under contracts no. DFG-Re856/6-l 
and DFG-Re856/6-2. 



For completeness, we write the explicit form of the various field equations of motion represented by Eq. (2.15): 



6. SUMMARY AND OUTLOOK 
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APPENDIX A: EXPLICIT FORM OF THE FIELD EQUATIONS OF MOTION 




P%Z[J] 




(A.l) 




(A.2) 



(A.3) 
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V a x Z[J] = J 2*&exp{tS + zS s } {V 2 x c a x - gf ahc V lx A\ x c%} . (A.4) 
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